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1. Introduction. 



r 25 1 

Deterministic Lanchester-type equations of warfare ( see TAYLOR and BROWN , 

[27] 

WEISS ) play an important role in military operations research for developing insights 

r 4 1 

into the dynamics of combat ( see , for example, BONDER and FARRELL , or BONDER and 

HONIG^), even though combat between two opposing military forces is a far more complex 

[9] 

random process. The classic Lanchester theory of combat ( see DOLANSKY ) considered 
constant attrition-rate coefficients. New operations research techniques for fore- 
casting temporal variations in fire effectiveness (caused by, for example, changes in 
force separation, combatant postures, target acquisition rates, firing rates, etc.) 
have generated interest in variable-coefficient combat formulations. Unfortunately, 
the resultant differential equations are not well studied. 

In this paper we present a mathematical theory for predicting battle outcome 
from initial conditions without explicitly computing force-level trajectories for variable 
coefficient Lanchester-type equations of modern warfare for combat between two homogeneous 
forces ( see Note 1) . The determination of conditions on initial values that predict 
(in the sense of necessary and/or sufficient conditions) force annihilation in such 
Lanchester-type combat leads to some new mathematical problems in the theory of ordinary 
differential equations. This force annihilation problem may be viewed as either a 
problem of determining the asymptotic behavior of the solution (depending on given 

initial conditions) or a problem of determining the range of the quotient of two linearly 

r 25 ] 

independent solutions to, for example, the X force-level equation 1 ( see Note 2). 

In either case, the classic ordinary differential equation theories ( see , for example, 
HILLE^^, INCE^*^, OLVER^^) are inadequate to supply all the answers sought. We 
show that questions of force annihilation can be reduced to the study of certain "expo- 
nential-like” Lanchester functions and may be simply answered by examining certain 
inequalities involving the initial conditions and by possibly consulting tabulations 
of new special functions that are suggested here. Our general results apply to a wide 
class of attrition-rate coefficients (namely, those that yield continuous force-level 
trajectories) . 



Thus, in this paper we provide a general theoretical framework for determining 
force annihilation without explicitly computing force-level trajectories for variable- 
coefficient Lanches ter- type equations of modern warfare. [Other modes of battle ter- 
mination are briefly discussed.] We introduce a canonical auxiliary parity-condition 
problem for such determinations. New exponential-like general Lanchester functions 
arise from the solution to this problem, and tabulations of these would facilitate 
force-annihilation prediction. Different mathematical forms for attrition-rate coeffi- 
cients lead to different auxiliary parity-condition problems. Our theory is applied 
to general power attrition-rate coefficients: exact force-annihilation results are 

obtained for cases of f, no offset 11 (modelling, for example, weapon systems with the 
same maximum effective range) ; and although qualitative results are obtained, future 
computational work is required for quantitative results for cases of "offset" (modelling 
for example, weapon systems with different maximum effective ranges). 

2 . Lanches ter *s Classic Formulation . 

r 13 1 

F. W. LANCHESTER hypothesized in 1914 that combat between two military forces 
could be modelled by 

dx/dt = -ay, dy/dt = -bx, (1) 

with initial conditions 

x(t=0) = Xq , y(t=0) = y Q , (2) 

where t = 0 denotes the time at which the battle begins, x(t) and y(t) denote 

the numbers of X and Y at time t, and a and b are nonnegative constants which 

are today called Lanchester attrition-rate coefficients and represent each side's fire 

effectiveness. Lanchester ( see McCLOSKEY^^ for his influence on operations research) 

considered this model (1) in order to provide insight into the dynamics of combat under 

"modern conditions" and justify the principle of concentration ( see Note 3) . We will 

accordingly refer to (1) as Lanchester 's equations of modem warfare . Various sets of 

physical circumstances have been hypothesized to yield them: for example, (a) both 

[27] 

sides use aimed fire and target acquisition times are constant ( see Weiss )> or (b) 



2 



both sides use area fire and a constant density defense ( see BRACKNEY ) • Other forms 
of Lanches ter- type equations appear in the literature, but we will not consider these 
here ( see Dolansky^ and TAYLOrJ^). 

From (1) Lanches ter deduced his classic square law 

b(x 2 -x 2 (t)) = a(y 2 -y 2 (t)) . (3) 

Consider now a battle terminated by either force level reaching a given "breakpoint” 

BP 

( see Note 4): for example, Y wins when x f = x(t f ) = Xp p = f x x Q but y f > y Bp 

BP 

= f v y , where t-, y- denote final values and X-^ denotes X f s breakpoint 

l U I L L Dr 

BP 

which is given fraction f of his initial strength. It follows from (3) that 

A 

Y wins » x Q < [(l-(fY P ) 2 }/{l~(f^ P ) 2 )1 1 ^ 2 v'a/b y Q , (4) 

which for a f ight-to-the-f inish becomes the classic result 

Y wins f ight-to-the-f inish ** x^ < /a/b y^. (5) 

Since, unfortunately, no relationship similar to (3) holds in general for variable 
attrition- rate coefficients, we observe that (4) may also be obtained from the time 
history of the X force level 

x(t) = { (xQ-yQ/a/bJexpC/ab t) + (x 0 +y Q /a/b) exp(-</atT t)} /2 , (6) 

BP 

via determining the time for X to reach his breakpoint (i.e. x(t=t ) = x ) 

X Dr 

t PP = (1/ t^JilnCt-Xgp+i/x^p+y^a/b-x^/tyQt/a/b - x Q }) , (7) 

BP BP 

and requiring t x < ^Y * The ^ e Y resu ^ t f° r obtaining (7) is that one of the two 
linearly independence solutions to the X force-level equation d 2 x/dt 2 - abx = 0 is 
the reciprocal of the other. For a f ight-to-the-f inish (7) becomes 

= (1/ (2v / ab) Hn({y 0 ' / a/b + x Q }/{y 0 /a/b - x Q }) , (8) 

where t denotes the time to annihilate the X force. We observe that (5) is an 
immediate consequence of (8) . 



3 



In many applications ( see Section 3 below) , one is interested in whether the 



battle will be terminated within a given time t^ . In this case Xq < /a/b is a 



necessary condition for X to be annihilated and annihilation occurs when t £ t . 

x g 



Thus, determination of whether force annihilation will occur within a given time involve 
consulting a tabulation of a transcendental function, here the natural logarithm ( see 
also (10) below). Similar results hold for other fixed force-level breakpoints. 

The time history of the X force level may also be written as 



x(t) = Xq cosh v / a^b t - y^ / a / b sinh /ab t . 



(9) 



[25] 



Taylor and Brown take (9) as their point of departure for a mathematical theory for 



solving variable-coefficient formulations. (7) does not follow directly from (9), but 



(5) does via x(t=t ) = 0 and 

A 



= (1/ Sab) tanh 1 (/b/a x Q /y 0 ) , 



( 10 ) 



since the range of the hyperbolic tangent is [0,1] for nonnegative arguments. 

The purpose of this paper is to generalize the above to the general case of 

variable attrition-rate coefficients. We use (6) as our point of departure and base 

our development on the observation that (5) follows directly from (6) , since the second 

term in brackets is always positive and goes to zero as t->“+ 00 . [We will ignore the 

physical impossibility of negative force levels in developing results like (5).] Thus, 

by (6) lim x(t) = -°° « x^ < Sajb yQ , whence follows (5). 
t-> 00 



3. Variable Attrition-Rate Coefficients . 

The pioneering work of S. Bonder^’ ^ on methodology for the evaluation of 

military systems (in particular, mobile systems such as tanks) has generated interest 

in variable-coefficient Lanchester-type equations and has led to improved operations 

[21 [3] 

research techniques for the prediction of such coefficients ( see BONDER * ; back- 

r 25 ] 

ground and further references are given in Taylor and Brown ). Thus, we consider 

dx/dt - -a(t)y, dy/dt = -b(t)x, (11) 

where a(t) and b(t) denote time-dependent Lanchester attrition-rate coefficients. 



4 



These coefficients depend on such variables as firing doctrine, firing rate, rate of 
target acquisition, force separation, tactical posture of targets, etc, ( see reference 
4). Without loss of generality, we may take a(t) = k^g(t) and b(t) = k^h(t) , where 
g(t) and h(t) denote time-varying factors such that a(t)/b(t) = k^/k^ = constant 
for g(t) = h(t). We will also refer to (11) as the equations for a square-law attri - 
tions process , since an ’’instantaneous 11 square law holds even when a(t)/b(t) is not 
constant ( see TAYLOR and PARRY ; also references 21, 22, and 24), 

A large class of combat situations of interest can be modelled with the following 
attrition-rate coefficients ( see reference 4) 

a(t) = k a (t+C) W and b(t) = k b (t+C+A) V , (12) 

where A,C ^ 0. We will refer to these coefficients as general power attrition-rate 

[25] 

coefficients . The modelling roles of A and C are discussed in Taylor and Brown. 

We will refer to C as the starting parameter , since it allows us to model (with 
y,v ^ 0) battles which begin within the maximum effective ranges of the two systems. 

We will refer to A as the offset parameter , since it allows us to model (again, with 
y,v £ 0) battles between weapon systems with different effective ranges. For example, 

let us consider Bonder's ^ constant-speed attack on a static defensive position ( see 
also references 22 and 25). Then we have 

dx/dt = -a(r)y = -a Q (l-r/R a ) y y , dy/dt = -3(r)x = -3 Q (l-r/Rg) V x , (13) 

where y,v ^ 0 and R^ denotes the maximum effective range of Y's weapon system 
(i.e. a(r) = 0 for r > R^) , and these parameters are given by (we assume that 

*6 ^ R a 2 V 

A = (R -R )/v, and C = (R -R rt )/v, (14) 

P Ot CL U 

where R^ denotes the battle's opening range and v > 0 denotes the constant attack 
speed. By considering (14) and Figure 1, the reader should have no trouble in under- 
standing our terminology for A and C. In this model (13) y, for example, is used 
to model the range dependence of Y's attrition-rate coefficient ( see Figure 2). 
Negative values of y and v are discussed in reference 25. Range is related to time 
by r(t) = Rg-vt. Since zero force separation is reached at time t^ = R^/v, the 
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maximum effective range of the weapon system and constant kill capability at zero range. 
[Notes: 1. The maximum effective range of the system is denoted R = 2000 meters. 

2. a(r=0) = a Q = 0.6X casualties/ (unit time x number of Y units) denotes the Y force 

weapon system kill rate at zero force separation (range). 3. The opening range of 
battle is denoted as R_ = 1250 meters and (as shown) R n < R . ] 



above model (13) provides interest in predicting battle termination within a given 
time t ( see Section 2 above) . 

Almost all previous work on the variable-coefficient equations (11) has developed 
infinite series solutions for force-level trajectories or represented these by tabulated 
functions ( see reference 25, in particular Section 3). Relatively little attention has 
been given to determining the qualitative behavior of solutions to (11) (such as predic- 
tion of battle outcome) without explicitly computing battle trajectories. Bonder and 
[4] 

Farrell ( see Note 5), however, have considered force annihilation within a given 

time. Using comparison techniques from the theory of ordinary differential equations 

r o 1 

( see , for example, CODDINGTON and LEVINSON 1 J ), they obtained a rather strong sufficient 
condition for the special case of the model (13) with y = v = 1 and A > 0. 



4. A Mathematical Theory for Conditions of Force Annihilation . 

Motivated by the constant-coefficient results ( see Note 6) , we introduce the 

, and E^, defined by 



exponential-like general Lanchester functions ' E^, E^, E^ 



dE^/dt = /k^/k^ a(t)E^ 
dE^/dt = /k~7k^ b(t)E^ 

and 

dE'/dt = -/k^TlT a(t)E^ 
dE^/dt = -/k~7k^ b(t)E x 



with E x (t=t Q ) = 1/Q, 
with E^(t=t Q ) = 1, 



(15) 



with E x (t=t Q ) = 1, 
with E^(t=t 0 ) = Q, 



(16) 



X Y X 

where t^ = max(t^,tQ), and t^ denotes the largest finite singular point ( see Taylor 
and Brown^*^; also p. 69 of Ince^^) ( see Note 7) on the t-axis for the X force- 
level equation 



d 2 x/dt 2 - { [l/a(t) ]da/dt}dx/dt - a(t)b(t)x = 0. (17) 

By assumption, then, a(t) and b(t) are positive continuous functions V t > t^. 

Both E^(t) and E^(t) satisfy (17). Since they are a fundamental system of 
solutions, we may use them to construct all solutions to (17). Thus, the solution to 
to (17) which satisfies the initial conditions (2) is given by 
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( 18 ) 



x(t) = { [x Q E^(t=0) - A a /k b y Q E x (t=0) ]E^(t) 

+ [x 0 Ej(t=0) + y 0 E^(t=0)]E“(t)}/2, 

and from (11) we obtain 

y(t) = { [y Q E x (t=0) - *^7kJ x Q E^(t=0) ]E^(t) 

+ [y Q E x (t=0) + x 0 E x (t=0) ]E y (t) }/2, 

where we have made use of the easily verifiable fact that ( see reference 25) 

E x (t)E^(t) + E x (t)E^(t) = 2 V t. 

Without further restrictions on a(t) and b(t), the systems (15) and (16) 
might not have a solution: it might be that one or more of the general Lanchester 

functions is unbounded at t^ ( see Section 4 of reference 25). For example, unbounded- 
ness occurs for power attrition-rate coefficients with p + v+2 = 0 and A = 0 ( see 

f 25 1 

equation (37) of Taylor and Brown ). Thus, we introduce 

r t 



(19) 



( 20 ) 



CONDITION (A) : 



tk V 



r 



a(s)ds and | b(s)ds are bounded for all finite 

t. 



Then, we have 



THEOREM 1: Condition (A) is a necessary and sufficient condition for (15) 



and (16) to have a continuous solution for all finite t^t^. 



PROOF: Sufficiency follows by Theorem 6.4.2 on p. 226 of Hille 



[ 10 ] 



(see also p. 64 



of LEE and MARKUS ^"^). To prove necessity, it suffices to consider E^~(t) . If a 

A 



continuous solution exists, then from (15) we have dE^/dt ^ /k^/k^ a(t), whence 

' t 

a(s)ds, and the theorem follows. Q .E .D . Let us henceforth 
t^ 



Ey(t) fel/Q+ /k, /k 
X b a 



assume that Condition (A) holds. 

Let us now consider how to choose the general Lanchester functions I and E^ 



so that they play the roles of a decaying exponential and an increasing exponential, 
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respectively, in the solution (18) to the variable-coefficient X force-level equation 

(17) [whence our notation of E and E^j • We recall the constant-coefficient result 

A A 

(6) and its consequence (5) obtained using lim exp(-/at7 t) = 0 and exp(-/abT t) > 0. 

t-H-oo 

For any positive value of Q it is clear that E+(t) , which satisfies (15) , grows 

A 

without bound just as an increasing exponential does. We will now show how to choose 

E so that it corresponds to a decaying exponential. Similar statements hold for 
A 

E^ and 

positive for all t so that by (16) they continuously decrease. Furthermore, we will 
be able to specify such behavior for E^ and E^ by our selection of the parameter 
Q in the initial conditions for (16) . 

The solution E (t), E (t) to (16) depends continuously on the parameter Q 
A x 

of the initial conditions ( see , for example, Hille^^). We denote this dependence 
by E x (t;Q), E^(t;Q). Let Q* = Q*(a( t) ,b ( t) ) denote the value of Q such that 

E x ( t ;Q=Q*) , E^(t;Q=Q*) > 0 for all finite t ^ t^. (21 

It follows from (16) that 

lim E“(t;Q*) = lim E~(t;Q*) = 0. (22 

It is intuitively obvious that such a Q* exists, and we will prove its existence 
in some particular cases below. As we shall see, knowledge of Q* provides valuable 
information about the qualitative behavior of force-level trajectories for the Lan- 
chester-type equations (11) . Let us refer to the problem of determining Q* such 
that (22) holds as the auxiliary parity-condition problem . Unless explicitly stated 
otherwise, for convenience we will denote, for example, E x (t;Q*) as E+(t) . 

Comment 1 : For a constant ratio of attrition-rate coefficients, i.e- 

a(t) = k^h(t), and b(t) = k^h(t) , (23 

where h(t) denotes the common time-varying factor of the two coefficients, it readily 
follows from the results given in references 4 and 20 that Q* = 1 and 



E^. Considering (16), we see that we should choose E^ and E^ to remain 
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and 



(24) 



E^(t) = E^(t) = exp{<J/(t) ) » 



E x (t) = ( t) = exp{-^(t)>. 



where i^(t) = /k^k^ J h(s)ds. 



Comment 2 : By Theorem 1 of Taylor and Brown 

form of (6) only if (23) holds. 



[25] 



, the solution (18) simplifies to the 



The above exponential-like general Lanchester functions may be related to Taylor 



[25] 

and Brown’s hyperbolic-like general Lanchester functions. Let 



C x (t) = x^t), s x ( t ) = x 2 ( t ^» Cy(t) = y 1 (t) , and S y (t) = y 2 (t), (25) 

where x^, x^, y^, and y^ denote the general Lanchester functions introduced by 
Taylor and Brown. Then, similar to the well-known relationships between the hyperbolic 
and exponential functions, we have 



C x (t) = (Q*E x (t) + E~(t)}/2, 



S x (t) = (Q*E x (t) -E”(t)}/(2Q*), 



(26) 



CyCt) = (Q*E^(t) +E y (t)}/(2Q*), S Y (t) = {Q*Ej(t) -E Y (t)}/2. 

The determination of Q* will be slightly simplified for general power attrition- 
rate coefficients (12) by considering a modified auxiliary parity-condition problem. 

For this purpose we introduce the new independent variable 

rt 



- K/ V i 



a(o)da , 



(27) 



and define s^ = s(t=0)^0 for t^^O. K is an, at present, undetermined parameter. 



It will be chosen so that a more convenient canonical system of differential equations 
arises in the modified auxiliary parity-condition problem. By Condition (A) , the 



transformation (27) is well defined for t^t^. It has an inverse t(s) , since a(t) > 



V t > t n . Letting 



4<s) =KE^(t(s)), e*(s) = E^(t(s)), e x ( s ) = E x (t(s)), and e Y (s) = E Y (t(s))/K, 



the substitution (27) transforms (15) through (17) into 



de x /ds = e Y 



with e x (s=0) = 1/Z , 



de^/ds = I(s)e x 



with 



e y (s=0) = 1, 



(29) 
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de x /ds = -e Y 
de^/ds = -I(s)e x 

and 

d 2 x/ds 2 - I (s)x = 0, 



where for any Q 



Z = Q/K, 



with e (s=0) = 1, 
A 

with ( s = 0) = Z, 



(3o; 



(31; 



02: 



and 



I(s) = ({b(t)/k fe }/{a(t)/k a })/K 2 , 



( 33 ; 



[ 12 ] 

is the invariant of the normal form (31) ( see p. 119 of KAMKE ) and 
(27), The parameter K will be chosen to simplify the form of I(s) . 
work the equation (31) will be easier to analyze than (17) • 

We will refer to the problem of determining Z* = Z*(a( t) ,b ( t) ) 



t = t(s) by 
In our later 

such that 



e (s;Z=Z*), e (s;Z=Z*) > 0 for all finite s ^ 0, 

A i 



( 34 ; 



as the modified auxiliary parity-condition problem . By (32) we then have that Q* = KZ : 
We also observe that the solution to (31) which satisfies (2) may be expressed 
in terras of these exponential-like general Lanchester functions for s ^ s^ as 

x(s) = { [x Q Ke Y (s=s 0 ) - >/k^7k^ y 0 e x (s=s 0 ) ]e x (s)/K 

+ [x Q Ke^(s=s 0 ) + y 0 e x (s=s 0 ) ]e x (s)/K}/2, (35; 

and similarly 

y(s) = {[y Q e x (s=s 0 ) - A b /k a x 0 Ke Y (s=s Q ) ]e^(s) 

+ [y Q e x (s=s 0 ) + /k^TiT - x Q Ke Y (s=s 0 ) ]e Y (s)}/2, (36] 

where by (20) and (28) 

e x (s)e Y (s) 4- e x (s)e^(s) = 2 Vs. (37) 

From our choice of Q* such that (21) and (22) hold, we can immediately infer 
the behavior of the solution (18) to (11) as t +® and similarly for y(t) (19). 



12. 



Thus, we have 



THEOREM 2: lim x(t) = -» if and only if 

x 0 E”(t=0;Q*) < A a /k^ y 0 E”(t=0;Q*). 

Equivalently, we may state 

THEOREM 3: Consider combat between two homogeneous forces described by 

(11). Assume that (11) applies for all time and that Y "wins" when 
x(t f ) = 0 with y(t f ) > 0. Then, Y wins if and only if XgE^t^OjQ*) < 

■'V \ y o E x (I;-05Q * > - 

Thus, we see that tabulations of such new exponential-like general Lanchester functions 

E x (t;Q*) and E^(t;Q*) would facilitate force-annihilation prediction. By our choice 

of tg only one such tabulation is necessary for given attrition-rate coefficients 

a(t) and b(t) ( see Note 7). Alternatively, we may express the force-annihilation 

r 25 1 

condition of Theorem 3 in terras of Taylor and Brown* s hyperbolic-like general 
Lanchester functions (see equations (25) and (26)). We have then 

THEOREM 3*: Consider combat between two homogeneous forces described by 

(11). Assume that (11) applies for all time and that Y "wins" when 
x(tf) = 0 with y(t^) > 0. Then, Y wins if and only if x^tQ*^^^) - 
S Y (t=0)} < ✓k a /k b y Q {C x (t=0) -Q*S x (t=0)}. 

As an immediate corollary to Theorem 3 we have 

COROLLARY 3.1: For t^ = 0 , Y wins if and only if x qQ* < ^k^ / k^ y^. 

One shortcoming of our above development is that Theorem 3 and its corollaries 

are basically existence theorems for the annihilation of one side by the other at some 

unknown future point in time. As the enemy initial force level decreases towards 

parity (i.e. equality holding in (38)), the time required to annihilate a force becomes 

larger and larger. There is, in fact, no limit to how large it may become. Moreover, 

[4] 

there is a large class of tactically significant battles ( see Bonder's constant- 
speed attack on a static defensive position discussed in Section 3 above) which has a 
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built-in time limit, denoted as t . Hence, it would be desirable to have a method for 

8 

determining from initial conditions (without explicitly computing the entire force-level 

trajectories) whether or not force annihilation will occur within a given finite time 

tg. For our general model (11), Theorem 3 tells us that 

x 0 E^(t=0;Q*) < /kVkh y Q E" (t=0;Q*), (38) 

is a necessary condition for x(t^) = 0 with t^ £ t , where tf. denotes the time at 

X X g X 

which the X force is annihilated. Motivated by the well-known constant-coefficient 
result (8) , one intuitively sees that determining whether or not t^ £ t will require 
the appropriate tabulations of new transcendents (i.e. new functions). Different such 
functions arise from different fundamental systems chosen to construct the solution to, 
for example, the X force-level equation (17). Let us now investigate which fundamenta 
system of solutions is the most useful. 

We begin our investigation by outlining for the new exponential-like general 
Lanchester functions how to determine from initial conditions (without explicitly com- 
puting the force-level trajectories) whether or not force annihilation will occur within 
a given finite time t^. Looking at (18) and setting x(t) = 0, we see that we must 
solve for E*(t)/E (t) = n(t) (see Note 8). For any other fundamental system (i.e. pai 
of solutions, we must still solve for such a quotient (cf. the constant-coefficient 
results (8) and (10)). Thus, our force-annihilation determination requires the use of 
tabulations of the quotient of two linearly independent solutions to, for example, the 
X force-level equation. Our question is now which quotient is the most useful (i.e. 
which fundamental system yields the most useful quotient) . We will show that the quo- 
tient of two exponential-like general Lanchester functions is not numerically satis- 
factory for such determinations, whereas the quotient of two hyperbolic-like general 
Lanchester functions is. 

Motivated by the result for a constant ratio of attrition-rate coefficients (23) 
that n(t) = exp{2ij>( t) } , we use the following notation for the n(t) of the X force- 
level equation 
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( 39 ) 



E 2X (t) = E x (t)/E X (t) ’ 

with being similarly defined. Assuming that (38) holds, we see from (18) 

that if x(t=t^) = 0, then 

E 2 X ( tx> = y 0 E ^(t=0;Q*) +x 0 E^(t=0;Q*)}/ 

y 0 E~(t=0;Q*) - x Q E~( t=0 ;Q*) } . 

We will prove below that E^CO = a str ^- ct ^y increasing function of t 

with initial value 1/Q* at t = t^. Consequently, the inverse function E^ ^(O = 
L^xCS) is a well-defined function of £V £ e [l/Q*,+°°). Hence, we may write that the 
time to annihilate X is given by ( see Note 9) 

C x = E 2x" 1({/k a /k b y 0 E x (t=0; Q* )+x 0 E Y (t=0;Q * )}/ 

y Q E^(t=0;Q*) -x 0 E^(t=0;Q*)}). (40) 



Considering the above developments, one may show that except when (23) holds, it is 
not possible to determine in the manner of (40) the battle-termination time when X’s 
breakpoint is positive, i.e. x^p > 0. 

We now show that E^CO (defined by (39)) is a strictly increasing function. 



We readily compute using (15), (16), and (20) that 
dE+ x /dt = 2/ylT a(t)/{E~(t) } 2 
whence follows the monotonicity. 



with E 2x (t=t Q ) = 1/Q*, 



(41) 



[Sometimes it is convenient to express results in terms of the modified exponen- 
tial-like general Lanchester functions e*(s) , e^s) , e*(s) , and ey(s) , which are 

based on the transformed time variable s defined by (27) . The above results are 
readily extended to this case. In particular, e 2x^ s ^ = e X^ S ^ 6 X^ S ^ a strictly 

increasing function of s for s ^ 0, e 2X^ S= ^ = K/Q*, and lim e 2x^ s ^ = 

s -H-o° 

Thus, (I) determination of Q* = Q*(a( t) ,b ( t) ) , and (II) tabulation of 
E^(t) and E^CO would allow one to determine whether or not force annihilation 
occurs in such battles with finite time limit without explicitly computing the entire 
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force-level trajectories ( see Figures 3, 4, and 5 of reference 25) . Unfortunately, 
there is a serious drawback to considering E^Ct) and E^yCt) : accurate tabulations 

are difficult (in fact, essentially impossible for large values of t) to generate, 
since both functions are basically increasing exponentials so that any error in their 
initial value 1/Q* (which, in general, can be only numerically approximately deter- 
mined) becomes tremendously magnified over time. 

We may develop numerically satisfactory functions for prediction of force annihi- 
lation within a given finite time, however, by considering the hyperbolic-like general 

f 25 1 

Lanchester functions of Taylor and Brown . Let us therefore define 



T X ( t) = S x (t)/C x (t). (42) 

The functions T (t) and T v (t) are analogous to the hyperbolic tangent, to which 
X I 

they reduce for a constant ratio of attrition-rate coefficients. Considering (25) and 

equation (16) of reference 25, we see that, for example, T v (t) does not depend on Q*, 

since S (t) and C (t) do not. Thus, T (t) and T v (t) are numerically suitable 
XX X Y 

for determining whether or not force annihilation will occur within a given finite time. 

[251 

Using the results given in Table I of Taylor and Brown , we readily compute 



that 

dT x /dt = /k^/k^ a(t) /{C x (t) } 2 with T x (t=t Q ) = 0. (43) 

Hence, T (t) is a strictly increasing function, and its inverse T ^ is well defined. 
X X 

Let us now establish an upper bound for T (t). By (26) and (42), we have 

X 

T x (t) = (l/Q*){Q*E x (t) -E"(t)}/{Q*E x (t) +E~(t)}, 

whence it follows that for t ^ t^ 

0 £ T v (t) < 1/Q* , with lim T v (t) = 1/Q*. (44) 

t-H-oo 

Thus, our current investigation has yielded important information about the asymptotic 
behavior of hyperbolic-like general Lanchester functions. 
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3 . 3 

To determine t such that x(t ) = 0, we write the solution to (17) which 

A A 

r 25 1 

satisfies the initial conditions (2) as 1 J x(t) = x Q^Cy( t= 0) C^(t) - S^C t=0)S x (t) } 

- y^/k^/k^ (C x (t=0)S^(t) - S^(t=0) C^( t) } and find that when (38) holds, t^ is given by 

t x = T x 1({x o C Y (t=0)+y O /k a /k b S x (t=0)}/ 

{ y o\/kb C x (t=0) + x Q S Y (t=0) }) . (45) 

We observe that by (26) the argument of the inverse function T ^ in (45) belongs to 
the range of ( see (44)) when (38) holds ( see Theorem 3 f ). For t^ = 0, (45) simpli- 

fies to t^ = T^ (/k^/k^ x q/Yq) • Whether or not force annihilation occurs within a 
given finite time t^ then depends on whether or not t^ £ t . Thus, (I) determination 
of Q*, and (II) tabulation of T^(t) and T^(t) would allow one to determine (with- 
out explicitly computing the force-level trajectories) the time at which a side is 
annihilated. 



5 . Application to Power Attrition-Rate Coefficients . 

Let us now apply the above general theory to (11) with the general power attrition- 
rate coefficients (12). We observe that in this case t^ = -C, where C ^ 0. In order 
that Condition (A) holds we must have y,v > -1. 

As we have seen in Section 4, our theory of f orce-annihilation prediction depends 
on knowing Q*, the solution to the auxiliary parity-condition problem. For the power 
attrition-rate coefficients (12), it is more convenient, however, to determine Q* via 
the modified auxiliary parity-condition problem (30) ( see also (34)). Hence, we apply 
the transformations (27) and (28) to (16). For the coefficients (12), equations (30) 
take the form 



de x /ds = -e Y with e x (s=0) = 1, 

de^/ds = -s^(l-by/s a ) V e x with e^s^) = Z, 

where the parameter K in (27) is given by K = ( /k^k^,/ (y+l) ) ^ ; and we have 

p = (y+1) / (y-Fv+2) , q = 1-p, a = l/(y+l), 6 = (v-y)/(y+l), and 



( 46 ) 
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Y = A(v / 1T^/(p+l)) 2 ^ lJ+V+2 ^ . For p,v > -1, we have 0 < p, q < 1. 

After we have solved the above modified auxiliary parity-condition problem (i.e. 

determined Z = Z* for (46) such that (34) holds) , we have all the information required 

to determine without explicitly computing the entire force-level trajectories whether 

or not force annihilation occurs in battles modelled with (12), We may apply Theorem 

3 via (28) (possibly using (26)) to see who can be annihilated and use results such as 

(40) (or equivalent ones for hyperbolic-like power Lanchester functions ( see Taylor and 
[251 

Brown )) to see if force annihilation occurs in a given finite time (for example, 
for battles modelled by (13)). When C = 0 for the coefficients (12) (e.g. for the 
model (13), = R ^ R 0 from (14)), we have by Corollary 3.1 

U ot p 

COROLLARY 3.2: For combat between two homogeneous forces modelled by (11) 

and (12) with C = 0, the X force can be annihilated if and only if 
Z*Xq < (^ a k^/ (y+1) ) ^ ^A^/k^ y^, where Z* = Z*(Y>y> v ) is such that 
(34) holds for (46) . 

We will next give exact analytic results for cases of no offset (i.e. 

A = 0 =* Y - 0) and discuss the difficulties of determining Z* when there is offset 
(i.e. A,y > 0). Moreover, results for the special case of no offset help provide 
a lower bound for Z* in the general case. 

6. Results for Power Attrition-Rate Coe fficients with No Offset. 

When the offset parameter A = 0, 

de x /ds = -e^ 

de^/ds = ~ s&e x 

Solving (47) by successive approximations, we obtain 

p2k s k(e+2) /{ ^ j(j_ p ) } 

- z l p 2k s k(e+2)+1 /{ £ j(j+p)>, (48) 

k=0 j=l 



e (s;Z) = l 
X k=0 



equations (46) become for y,v > -1 



with e (s=0) = 1, 
X 



with e Y (s=0) = Z. 



(47) 
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I 



jtfhich is not the most useful result for large s. Without explicitly computing e^(s;Z) 

Lt is impossible to determine from (48) how e (s;Z) behaves for increasing s. However, 



let us write e^(s;Z) as [observe that e^ satisfies the generalized Airy equation 
( see SWANSON and HEADLEY 1 ) d 2 e^ / ds 2 - s^e^. = 0, which is well known to be reducible 
to Bessel’s equation] 

e~(s ;Z) = r(q)p" q {A (s) + [l-Zp q_P r(p)/r(q)] [ps 1/2 I (S) > , (49) 

A P P 

for any arbitrary Z, where A (s) denotes the generalized Airy function of the first 

p 

[19] 

kind of order 8 ( see Swanson and Headley ) , I (S) denotes the modified Bessel 

P 

( g+2) / 2 

function of the first kind of order p, and S = 2ps v . Also, 

e^(s;Z) = r(q)p' q {(2p/ir)(sinp7r) s (8+D/2 (g) - [l-Z p q ' P r(p) /r (q) ]df /ds(s) } , (50) 

where K (S) denotes the modified Bessel function of the third kind (also called 

1/ 2 

Macdonald’s function) of order q and F(s) = ps I (S) . 

[19] 

The behavior of A (s) for s ^ 0 is readily seen from (see Swanson and Headley , 

p 

A 0 (s) = ( 2p/7r) (sin pu) s^^K (S) . It is readily seen (see p. 119 and p. 123 of LEBEDEV^^ 

B p 

or pp. 250-251 of Olver^^) that K^(S) is strictly decreasing and positive, and 

lim s^K (S) = 0, where v is any real number. It is then clear ( see also p. 1404 of 
s-H~ P 

reference 19) that A (s) , I (S) > 0 V s ^ 0, lim A.(s) = 0, and lim I (S) = +°° . 

B P s-H~ 6 _ _ S -H~ P 

Consequently, the requirement that (34) holds (i.e. that e^ and e^ behave like 
strictly decaying exponentials) means that the second term in the expressions (49) and 
(50) must vanish. Thus, for y = 0 we have 



z*( Y =o,p,v) = P P q r(q)/r( P ) . 



(51) 



Hence, Theorem 3 becomes ( see Note 10) 

THEOREM 4: Consider combat between two homogeneous forces described by (11) 

with attrition-rate coefficients (12) with A = 0. Assume that the model 
applies for all time and that Y ’’wins” when x(t^) = 0 with y(t^) > 0. 

Then Y wins if and only if x^e^ (s=s^ ;Z*) < A^/k^ yg( A^k^/ (p+1) )^ P e^(s=s^ ; Z*) , 
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where Z* = Z*(y=0,y,v) = p P q r(q)/r(p), e x (s;Z*) = p P r(q)A^(s) , and 
e Y (s;Z*) = ( 2p / tt) (sin pir) p q T (q) s^ +1 ^\^(S) . For C = 0, we have = 0 
so that Y wins if and only if x^r(q) < /k^/k^ y qT ( p) ( /k^k^ / ( y+v+2) ) q P . 

We observe that the infinite series form (48) is not of any value for determining 

asymptotic properties of the solution e to (47) (and consequently Z*) , although it 

A 

is useful for computational purposes. 



7. Offset Linear Attrition-Rate Coefficients . 

When the offset parameter A > 0, explicit analytic results for Z* are 

apparently not possible. Before considering the general case of y,v > -1 and y > 0, 

we will find it instructive to consider the special case of offset linear attrition- 

T41 

rate coefficients (i.e. y > 0 with y = v = 1) studied by Bonder and Farrell . This 
examination will show us why analytic results for Z* are elusive in cases of positive 
offset. 

For y > 0 with y = v =1, equations (30) and (31) become 

d 2 x/ds 2 - (l+y/t/s)x = (52) 



with initial conditions x(s=0) = x^ and dx/ds(s=0) = -/k^/k^ y^ , and 
de x /ds = -e^ with e x (s=0) = 1, 

de^/ds = -(l+y/v // s)e x with e^(s=0) = Z. 



(53) 



Let us now consider solving the above modified auxiliary parity-condition problem 
(i.e. determining Z = Z* for (53) such that (34) holds). Solving (53) by successive 
approximations, we may write 



e x (s) = h(s) - Zw(s) , and ^(s) = ZH(s)-W(s), (54) 

where h(s) = h(s;y) = h(A=s , y=y/v's~) , w(s) , H(s) , and W(s) denote the auxiliary 

r 25 ] 

offset linear Lanchester functions introduced by Taylor and Brown 1 J . Infinite series 
representations of these functions are given in reference 25. Subsequent research has 
shown that these hyperbolic-like Lanchester functions possess the properties given in 
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Table I. Unfortunately, information on the asymptotic behavior of the auxiliary offset 
linear Lanchester functions for large s > 0, which is needed to solve the modified 
auxiliary parity-condition problem (53) , is lacking and apparently not obtainable by 
the standard methods involving integral representation ( see Ince^^ or Olver^^). 
Consequently, we have not been able to develop an explicit analytic expression for 
Z*(y>0,p=l,v=l) , although we give upper and lower bounds for Z* in the next section 
for general y,v > -1. Additionally, we should point out that there are computational 
difficulties in searching for Z* via its definition (34): (a) one doesn’t know how 

large to take s for ’’satisfactory" results, and (b) numerical difficulties in evaluat- 
ing e^(s;Z) and e^(s;Z) as given by (54) occur for large values of s (since we 
are taking the difference of two very large numbers and, at least on a digital computer, 
can retain only a limited number of significant digits in these numbers) . 

Equation (52) is deceptively simple looking. Using variation of parameters, we 
may also express its solution as 



x(s) 



Xq cosh s - y^A^/k^ sinh s + y 



'S 

sinh (s-a)x(a) / Ado . 

0 



(55) 



Although (55) is a simple looking expression, this Volterra integral equation is, unfor- 
tunately, no easier to solve than (52) and leads to the same results as given by 
[22] r 251 

Taylor and Taylor and Brown . However, in the next section we show an easy 

way to obtain valuable information about Z* directly from (53) . 



8. Bounds on Z* for Power Attrition-Rate Coefficients With Positive Offset . 

We will now develop upper and lower bounds for Z*(y>0,y,v). These bounds 
establish the existence of Z* (and consequently Q*) for general power attrition- 
rate coefficients (12) by the continuous dependence of solutions to (46) on the initial 
conditions . 

The following two lemmas will be used to obtain an upper bound for Z*(y,p,v) 
for y , v > -1. 
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TABLE I. Properties of the Auxiliary Offset 



Note: 



Linear Lanchester Functions 



1. dh/ds = W, 



dW/ds = (l+y/v^h 



2. dw/ds = H, 



dH/ds = (1 -Py//s")w 



3. h(s)H(s) - w(s) W(s) =1 Vs 



4 . 



h(s=0) = H(s=0) = 1 



5 . 



w(s=0) = W(s=0) = 0 



6 . 



h(s;y=0) = H(s;y=0) = cosh s 



w(s;y=0) = W(s;y=0) = sinh s 



We use the notation h(s) = h(s;y) = h(A=s ,y=y / /s) 

r 25 1 

( see Taylor and Brown 1 J ), and similarly for w(s) , H(s) , and W(s) . 
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LEMMA 1: For 6^1 and x,y > 0, 2^ ^{x^+y^} ^ (x+y)^. 



PROOF: For 6 £ 1, f(x) = x is a convex function* A well-known theorem for convex 

functions says that { f (x) + f (y) } /2 ^ f([x-Fy]/2), whence follows the lemma. Q *E *D . 

p p p 

LEMMA 2: For 5^1 and x,y > 0, x+y ^ (x+y) . 

6 p x 

PROOF : Dividing by (x+y) , we need to show that [x/(x+y)] + [y / (x+y) ] ^ 1. If x 

and y > 0, then x/ (x+y) and y/(x+y) are < 1. Hence, for any 6^1 we have 

p p p 

[x/ (x+y) ] ^ x/(x+y) so that [x/(x+y)] +[y/(x+y)] ^ [x/(x+y) ] + [y / (x+y) ] = 1. Q.E.D . 

Using the above lemmas, we now prove Theorem 5. The upper bound to be given 
in Theorem 5 might be improved upon, although we feel that for computation determina- 
tion of Z* by interval search, it is not essential to have a better bound. 

THEOREM 5: For p > -1 and v 1 , we have Z*(y,p,v) < 1 + 2 V “ 1 (p+1) 2 / 

[ ( v+1) (p+v+2) ] +-y V 2 V 1 (p+1) 2 / (p+2) ; while for p > -1 and -1 < v £ 1, 
we have Z*(y,p,v) < 1 + (p+1) 2 / [ (v+1) (p+v+2) ] + y V (p+1) 2 / (p+2) . 

PROOF: Recalling (34), we have from the first equation of (46) that e^(s;Z*) < 1 

for s > 0; and from the second equation of (46) we then obtain for s > 0 and v ^ 1 

de^/ds > -s^(l+y/s a ) V -2 V ^s^ ( l+y V /s aV ) , (36) 

the latter inequality being a consequence of Lemma 1. From (30) we obtain 

e^(s) > z*- 2 v_1 s (v+1)/(p+1) • (p+1) /(v+1) - 2 V-1 y V (u+l)s 1/(lJ+1) . (57) 

Using (57) and considering the first equation of (46), we obtain e^(s;Z*) < U(s;Z*)> 
where 

U(s;Z*) = l+2 v_ h(p+l) 2 /(p+2)}{s (p+V+2)/(p+1) • (p+2) /[ (v+1) (p+v+2) ] 

+ y V S (P+2 ) /(p+ 1) } - Z *s . 

Since we must have 0 < e (s;Z*) < U(s;Z*) for all s ^ 0, it follows that for s = 1 

A 

we must have U(s=i;Z*) > 0, whence follows the theorem for v ^ 1. Lemma 2 and 
similar arguments are used to prove the theorem for -1 < v £ 1. Q • 
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Let us now consider the development of a lower bound for Z*(y>0,y,v). Before 
proving the key lemma (Lemma 3) for the proof of Theorem 6, we discuss some preliminary 
considerations . 

roz: i 

As shown by Taylor and Parry 1 , force annihilation for square-law attrition 
processes sometimes may be predicted by considering the force-ratio equation. For 
equations (46) , the "force ratio" u = e v /e satisfies the Riccati equation 

A i 

du/ds = s e (l-Pr/s a ) V u 2 - 1, (58) 

with initial condition u(s=0) = 1/Z. We observe that u(s^) = 0 if X is annihilated 
(i.e. e x^ s f^ = 0 but e^(s^) > 0), and u(s^) = if ^ is annihilated. For 
Z = Z* = Z*(y,y,v) such that (34) holds, we have a "draw," and u*(s) = u(s;Z=Z*) > 0 
and finite for all finite s ^ 0. Since (58) has a singularity at s = 0, we apply 
the transformation t = (s a +y) P+ l to obtain for x ^ Xq = y P+ ^ 

du/dx = x^u 2 - (l-y/x a ) P , (59) 

with initial condition u(x=y P+ ^) = 1/Z. 

We now develop a lower bound for Z*(y>0,y,v) by comparing results for y > 0 
with those for y = 0. For y = 0, we denote u as w and obtain for x ^ 0 

dw/dx = x^w 2 -l. (60) 

Corresponding to Z = Z* = Z*(y=0,y,v), we have via (49) and (50) 

w*(t) = T _e/2 K (T)/K (T) > 0, (61) 

P q 

( g+2) / 2 

where T = 2px v . Since K^(x) is finite and > 0 for all v,x > 0 ( see 

[141 

Lebedev , p. 136), we readily verify that w*(x) is finite for all x > 0 and 
that w*(x=0) = p^ P r(p)/T(q). Let us observe that w(x^) > w*(x^) = w(x^;Z*) 25 we 
have w(x) > w*(x) V x ^ 0 and W ( T 2 ^ = +°° f° r some finite T 2 > T l 9 si nce D = w - w* 
satisfies dD/dx = x P (w+w*)D. Consequently, w(x;Z) corresponds to Z < Z* and 
w(x;Z) becomes infinite at some finite time ( see equations (49) and (50)). We now 
state and prove the key lemma for developing a lower bound for Z* 
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LEMMA 3: Let w*(t) be given by (61) and let u(x) = u(t;Z) satisfy 
(59) for t ^ Tq = Then if u(t^) ^ w*(t^) , it follows that 

u(x) >0 V t £ and u(t£) = +" for some finite Conse- 

quently, Z*(y>0,y,v) > Z. 



PROOF: Consider D = u-w. It satisfies for T ^ t q v i- a (59) and (60) the equation 
dD/dr = t^(u-Fw)D+ {1-(1-y/t W ) ^} . If D(t^) = u(t^)-w(t^) ^ 0, then dDVdT(i^) > 0 
and D(t) > 0 Vt > T^, Thus, when u(t^) £ w*(t ) , we can find w(t^) for ^ 

such that u(x^) > w(t^) > whence follows the lemma from u(t) > w(x;Z) 

V t ^ t 2 with Z < Z*(y=0,m,v) and the above observation that w(t;Z) becomes 
infinite. Q.E.D . 

Letting x^ = x^, we °^ )t:a: ^ n 



THEOREM 6: Z*(y>0,p,v) > l/w*(y P+1 ) = 2 K (T_) /K (T_) , where 

U q U p U 



T o ■ 1 



-*- 1 and T 0 - 



Since w*(t=0) ^ w*(t) for g £ 0, we have as an immediate corollary 



COROLLARY 6.1: For v £ p, Z*(y>0,p,v) > p P q r(q)/r(p). 



Let us observe that for 6 > 0, the lower bound given in Corollary 6.1 is weaker than 
that in Theorem 6. 



9 . Future Computational Work . 

As we have seen above in Section 4, force-annihilation prediction depends on 
knowing the parity-condition parameter Q*, which may be called ,f the Y equivalent 
of an X force of unit strength.” We have explicitly determined Q* (via determining 
Z* for the modified auxiliary parity-condition problem (46)) for the power attrition- 
rate coefficients (12) in the case of no offset, i.e. A = 0. Tabulations of, for 
example, the new modified exponential-like general Lanchester functions e (s;Z*) and 

A 

e^(s;Z*) would facilitate force-annihilation prediction (see Theorem 4). It remains 
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to determine Q* for cases of positive offset, i.e. A > 0. As discussed in Section 7, 

analytic results for Z* in the modified auxiliary parity-condition problem (46) with 

y = A* ( /k^k^/ (y+1) ) > q are apparently not possible by the usual analytical 

methods so we must turn to numerical methods. It appears that a large number of cases 

[25] 

of tactical interest ( see Taylor and Brown ) would be covered by determining Z* 

for y,v = 0,1, 2, 3 and for a range of values of y > 0. One would be interested in, 

for example, plotting Z* versus y for fixed values of y and v. 

Since we have developed upper and lower bounds for Z* when y > 0 ( see 

Theorems 5 and 6) , we can use standard one-dimensional search techniques ( see , for 
r 28 1 

example, WILDE L J ) to calculate an approximate value of Z* with any predetermined 
degree of accuracy (of course, depending on how much computation we wish to do). Since 
(34) must hold for all s ^ 0, we must determine how long (i.e. for how large a value 

_ A _ A 

of s) to carry out computations of e^(s;Z) and :^(s;Z) in the modified auxiliary 
parity-condition problem (46) for a given trial value of Z* (denoted as Z) to see 
whether it is too large or too small. Although - e^ < 0 or e^ < 0 clearly terminates 
calculations for a given Z, it would be desirable to be able to terminate computations 
before this occurs, especially when Z is near Z*. In the future we will show that 
by considering the Ricatti equation (59) one can "cut off" conputations for a given 
value of Z well before either of the two annihilation conditions is actually reached. 
This acceleration in the termination of trial computations is particularly useful when 

A 

force annihilation occurs for large s (i.e. for Z near Z*) . Additionally, the 

qualitative behavior of solutions to (46) (or (11) , for that matter) is probably best 

understood by considering the Ricatti equation (59) (or, equivalently, (58)). 

As discussed at the end of Section 4, prediction of force annihilation within 

a given finite time involves the use of tabulations of the quotient of two linearly- 

independent general Lanchester functions. We have indicated in Section 4 that the 

hyperbolic-tangent-like Lanchester functions (e.g. T (t) as defined by (42)) are to 

A 

be preferred for reasons of the accuracy of their numerical computation. Thus, there 
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is a need for tabulations of T^(t) , whose range is [0,1/Q*) for te [0,-H»). For 

the power attrition-rate coefficients (12) with no offset, the power Lanchester functions 

(also called Lanchester-Clif f ord-Schlaf li (or LCS) functions), however, were inappro- 

f 251 

priately defined in Taylor and Brown to yield such tabulations. Thus, our newer 

theory of force-annihilation prediction, which also involves tabulations of canonical 

solutions (i.e. canonical Lanchester functions) to variable-coefficient Lanchester-type 

equations of modern warfare, has suggested some refinements in the definition of Taylor 
r 251 

and Brown’s auxiliary power Lanchester functions ( see Note 11). It would be desir- 
able then to redefine the LCS functions to fit within the framework of Section 4 and 
to develop tabulations of their quotient. If this were to be done, linear combat models 
with power attrition-rate coefficients (no offset) could be analyzed with somewhat the 
same ease as constant-coefficient linear models (i.e. (1)). 

10 . Extension to a More General Model . 

In this section we show that all our above force-annihilation results (except 
those for force-annihilation within a fixed finite time) also apply to a special case 
of a more general model. Moreover, comparison techniques may be used to extend these 
results in weakened form to the general case of this more comprehensive model. 

Let us consider the following Lanchester-type equations 



dx/dt = -a(t)y-g(t)x 


with x(t=0) = Xq, 


dy/dt = -b(t)x-a(t)y 


with y(t=0) = y Q , 



where a(t) , b(t), a(t) , and g(t) > 0. We may think of these equations as modelling, 

for example, aimed-fire combat between two homogeneous infantry forces with superimposed 
effects of supporting weapons (which are not subject to attrition and deliver area 
fire against the enemy infantry) ( see Taylor and Parry for a further discussion 

of this model (62)). In this case, a(t) and 3(t) are attrition-rate coefficients 

r 25 1 

which reflect the fire effectiveness of the supporting weapons 1 . Then, the force 
ratio u = x/y satisfies the generalized Riccati equation 
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du/dt = b(t)u 2 +{a(t) - f3(t)}u-a(t) 



(63) 



with u(t=0) = x q/y q • 

For equal effectivenesses of the supporting fires [i.e. a(t) = 0(t)], equation 
(63) simplifies to 

du/dt = b(t)u 2 -a(t), (64) 

which is the same Riccati equation satisifed by the force ratio for the model (11). 
Hence, when a(t) = 8(t) V t ^ 0, a battle’s outcome (in terras of the force ratio) 
is the same for the two models (11) and (62) , although the battle ends more quickly 
for (62). Thus, in this special case all our above results on force annihilation with- 
out time limitation (e.g. Theorems 2 through 4 and their corollaries) developed for (11) 
[in general or with the coefficients (12)] also apply to the more general model (62). 
Furthermore, comparison techniques ( see , for example, Hille^^) may be used to extend 
these results in weakened form to (62). Consequently, we see that the force-annihi- 
lation results developed in this paper are indeed of a fundamental nature. 

11 . Summary . 

We have presented a mathematical theory for predicting force annihilation for 

variable-coefficient Lanchester-type equations of "modern warfare" for combat between 

two homogeneous forces without explicitly computing force-level trajectories ( see 

[25] 

Note 12). This theory both generalizes and complements Taylor and Brown’s theory 

of canonical solution forms. Our force-annihilation theory provides guidance for 
certain parameter determinations and development of tabulations of Lanchester functions 
(beyond those suggested in reference 25) that would allow one to parametrically analyze 
variable-coefficient models with somewhat the same facility as constant-coefficient 
ones . 

We have shown that force annihilation can be predicted from initial conditions, 
without explicitly computing force-level trajectories, by knowing a parity-condition 
parameter Q* , which is the solution to a canonical auxiliary parity-condition problem. 
In general, this prediction would be facilitated by having tabulations of certain 



Lanchester functions available. In particular, to apply Theorem 3 it would be convenient 
to have tabulations of the new exponential-like general Lanchester functions E and 

A 

E^. However, if one wishes to invoke the equivalent Theorem 3 f , then tabulations of 
the hyperbolic-like general Lanchester functions C , S , C , and S would be 

A A I A 

desirable to have. The parity-condition parameter Q* was shown to be related to 

the range of the quotient of two hyperbolic-like general Lanchester functions introduced 

[25] 

by Taylor and Brown . Consequently, our force-annihilation theory not only provides 
new information about the mathematical properties of hyperbolic-like Lanchester functions 
but also provides guidance for selecting canonical Lanchester functions. 

We applied our general theory to the specific case of general power attrition- 
rate coefficients. Considering a modified auxiliary parity-condition problem, we 
explicitly determined Q* (via Z* of the modified problem) for power attrition-rate 
coefficients with no offset and gave upper and lower bounds for Z* for cases of 
positive offset. Consequently, in the future one -dimensional search techniques may 
be used to numerically determine approximate values of Z* as a function of the off- 
set parameter A. We finally showed that certain of our force-annihilation results 
also applied to a more general linear differential equation combat model. 

These results may be used in the analysis of the dynamic combat interactions 
between two homogeneous forces with time- (or range-) dependent weapon system capa- 
bilities. There is interest today in such analytic models because of improved opera- 
tions research techniques for predicting Lanchester attrition-rate coefficients, in 

particular their temporal variations ( see references 2 through 5 and 25) . Further 

[4] 

discussion of such applications may be found in Bonder and Farrell , Bonder and 
Honig^, Taylor , and Taylor and Brown^^. 
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NOTES 



1. Bonder and Honig^ point out, however, that force annihilation may not be the 



appropriate criterion for evaluating many military operations, expecially when force 

[4] 



annihilation does not occur. See pp. 192-242 of Bonder and Farrell J for a detailed 
Lanchester-type analysis of an attack scenario for which other "end of battle condi- 
tions'* play the major role in the evaluation process. Nevertheless, it is of interest 
(especially for developing insights into the dynamics of combat) to be able to easily 
predict (without computing force-level trajectories) the occurrence of force annihilation 
Such results are not only of intrinsic interest but also are useful in the optimization 



of combat dynamics, i.e. the combining of Lanchester-type models with generalized 



f2H F 23 1 

control theory (i.e. optimal control/differential game theory) ( see Taylor J,L J ) . 



r/1 TOOT r r 1 

2. Previous work by Bonder and Farrell 1 , Taylor 1 , and Taylor and Brown 1 shows 



that new transcendental functions arise even in the case of linear attrition-rate 
coefficients reflecting weapon systems with different effective ranges (i.e. the 



coefficients (12) with y = v = 1 and A > 0) . For example, the differential equation 

[12 



(52) could not be found among the 445 linear second order equations tabulated in Kamke 



Moreover, even when one can express a solution in terms of previously known transcendents 



the appropriate tabulations ( see , for example, ABRAMOWITZ and STEGUN ) may not exist 
( see Section 5 of reference 25) . As the equations of mathematical physics have pro- 
vided interest in many previously studied transcendents, so does the Lanchester theory 
of combat with time-varying fire effectivenesses provide interest in new transcendents. 



3. The influential 19th-century German military philosopher, Carl von Clausewitz 
(1780-1831), stated in his classic work On War ( Vom Kriege ) ( see p. 276 of reference 7), 
"The best Strategy is always to be very strong , first generally then at the decisive 
point. ...There is no more imperative and no simpler law for Strategy than to keep the 
forces concentrated." 
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4. As pointed out in reference 26, the entire topic of modelling battle termination 
is a problem area in contemporary defense planning studies, and there is far from 
universal agreement as to even which variables should be taken as the significant 

r 23 1 

variables for modelling this complex process. For further references see Taylor 

[4] 

5. Bonder and Farrell take range (i.e. force separation) to be the independent 

[22] [25] 

variable in their work, while Taylor and Taylor and Brown take time as we have 

done in this paper. 

6. Recalling the constant coefficient result (6), we consider 

d{exp(v / ab~ t) }/dt = a/b?a exp ( /ab" t) = b /a/b exp(/atT t) , 

and 

d{exp(-v / atT t) }/dt = -a/b/a exp(-v^TtT t) = -b/a/b exp(-/ab t) , 
to obtain motivation for (15) and (16) . 

7. We take a(t) and b(t) to be analytic in. the (finite) complex plane except for 

a finite number of singularities on the real axis. The singularities of (17) then 

occur at the zeros and singularities of a(t) and at the singularities of a(t)b(t). 

Consequently, t^ belongs to the set of points consisting of the zeros and singulari- 

[25] 

ties of a(t) and b(t) ( see Taylor and Brown ). We define t^ this way in order 
to reduce the number of tabulations of exponential general Lanchester functions required 
for force-annihilation analyses ( see Theorem 3). For example, for the general power 
attrition-rate coefficients (12) we have t^ = -C , and for fixed A ^ 0, p , and v 
only a single tabulation of e^. and e^ is required to handle all problems with 
C £ 0 ( see , for example. Theorem 4). 

8. It is well known ( see , for example, pp. 647-650 of Hille^^ or p. 120 of Karake^ ^) 

that the quotient of two linearly independent solutions to (31) , which is equivalent 

[18 ] 

to (17), satisfies Schwarz’s (third order) differential equation ( see SCHWARZ ) 
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(n ,s> = - 21 (s) , 



where q denotes the quotient of two linearly independent solutions to (31) [i.e. 

q = e*(s)/e (s)], { n * S } = q ,n - (3/2) (q") 2 /r| T denotes the Schwarzian derivative of 

X X 

q with respect to s, q 1 denotes dq/ds, etc., and I(s) denotes the invariant 
of the normal form (31). For numerical computation of q, however, there is no advan- 
tage to consider this third order equation, and it is preferable to calculate q from, 

for example, q(s) = e*(s)/e (s), where e+ and e also satisfy (29) and (30). 

XX XX 

9. Recalling our development of (7), we see that, except for the special case in which 
(23) holds, this same approach fails to yield the time for X to reach his breakpoint 
(assumed to be positive) [i.e. t^ such that x(t=t^ ) = x^p>0]. Consequently, it 
is apparently impossible to predict in the manner described in the main text the out- 
come of a fixed force-level breakpoint battle with positive breakpoints unless (23) 
holds . 



10. For the case of power attrition-rate coefficients with no offset (i.e. A = 0 in 

(12)), the second annihilation condition given in Theorem 4 (i.e. the one for C = 0) 

and an equivalent form of the first (i.e. the one for C ^ 0) may be developed by 

inspection when one expresses, for example, the time history of the X force level 

(which satisfies (17)) in terms of the so-called generalized Airy functions ( see Swanson 
U91 

and Headley ). For example, for C = 0 we have 

x(t) = (p P /(2p)){X+Y}A e (T) + (p P /(2v£)){X-Y}B g (T) , 



where A and B denote the generalized Airy functions of the first and second kinds 
p p 

of order g, X = x Q r(q) , Y = y Q A a /k b r(p) (/k^/ (p+v+2)) 1 2p , and T = 

( A^k^/ (q+1) ) 2p t P+ ^. The result given in Theorem 4 for C = 0 follows from the proper- 
ties of the generalized Airy functions (i.e. A (£) , B (£) > 0 V £ ^ 0, lim A (£) = 0> 

v v 



and lim B (?) = +“) and the above representation for x(t). Unfortunately, this result 
V 

does not generalize to other cases of interest, although it did motivate our general 
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theory of force annihilation developed in this paper. The generalized Airy functions 
may be considered to be generalizations of the exponential function ( see p. 446 of 
reference 1 or p. 393 of reference 17 for plots of the standard (i.e. 8 = 1) Airy 

functions) and arise in the study of the asymptotic behavior of solutions to certain 
differential equations. 

11. Although theoretically results for power attrition-rate coefficients with no 

offset are expressible in terms of "known" transcendental functions, new Lanchester 

[251 

functions were introduced by Taylor and Brown because of lack of tabulations of 
these in many cases of interest. 

[91 

12. In his well-known survey paper on the Lanchester theory of combat, Dolansky 
suggested the development of outcome-predicting relations without solving in detail 
as one of several problems for future research. 
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